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Abstract 

The spin-Peierls instability in spin 1/2 XY chain coupled to dispersionless phonons of frequency u) has been studied 
in the nonadiabatic limit. We have chosen the Lang-Firsov variational wave function for the phonon subsystem to 
obtain an effective spin Hamiltonian. The effective spin Hamiltonian is then solved in the framework of mean-field 
approximation. We observed a dimerized phase when g is less than a critical value and an anti-ferromagnetic phase 
when it is greater than a critical value . The variation of lattice distortion, dimerized order parameter and energy 
gap with spin phonon coupling parameter has also been investigated here. 

1. Introduction 

The discovery of the spin Peierls transition in CuGeO^ has sparked an intense effort to study the properties 
of this quasi one-dimensional magneto-elastic system where the coupling of the magnetic to the lattice degrees of 
freedom leads to a phase transition into a dimerized phase. This magneto elastic transition is due to the competition 
between the gain in magnetic energy due to dimerization and the loss in the elastic energy of the lattice distortion. 
Recently quite a large number of experimental [^]-^| and theoretical |6|-|l3| works have been performed to investigate 
the various aspects of the spin Peierls systems. Except a few 0,^5| most of the theoretical investigations rely on 
the adiabatic treatment of the phonons. In adiabatic approaches one assumes that the phonons responsible for the 
distortion have low energy with respect to the characteristic energies for the spin systems (e.g. the gap=A). The 
experimental evidence for CuGeO^ indicates that the application of the adiabatic approximation to these system is 
not sufficient. Regnault et al. investigated the spin dynamics of spin Peierls system CuGeO^ by inelastic neutron 
scattering. Their result confirmed the existence of a gap (A) in the magnetic excitations at A = 2meV = 23K and the 
Heisenberg exchange parameter was considered to be Ji — 10.6meV — 115A'. Braden et al. [fl7| found by symmetry 
that four optical phonons are possible candidates for the spin Peierls distortion in CuGeO^. Among these four modes 
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two, one of energy 330°K and other of energy 150 o -ftf are experimentally found to be the suitable candidates for the 
spin Peierls distortion. In both cases we find that the phonon frequency is larger than J as well as A. When u> > A 
the spin phonon interaction is unretarded and non-adiabatic effects or quantum lattice fluctuations become important. 
Fradkin and Hirsch |Q considered the Su-Schrieffer-Heeger model of electron-phonon interaction for spin-less 
fermions and spin ^ electrons in one dimension to investigate the stability of the Peierls-dimerized ground state against 
quantum fluctuations. In this work they have shown by renormalization group arguments and by quantum Monte 
Carlo simulation that for spin-less fermions quantum lattice fluctuations destroy the long-range dimerization order 
when the fermion-phonon coupling constant is small and predicted a transition from an undimerized ground state 
to a dimerized phase when the fermion-phonon interaction is larger than a critical value. Campbell and Bishop pQ] 
independently confirmed the findings of Fradkin and Hirsch. Recently Caron and Moukouri fHf] studied the XY spin 
chain coupled to dispersionless phonons by the density matrix renormalisation group (DMRG) method and showed 
that quantum fluctuations reduce the spin Peierls gap and even destroy the dimerization when the phonon frequency 
gets appreciably larger than the gap. 

The coupled spin-phonon system for all values of the coupling parameter is a very difficult problem. However, one 
may gain considerable insight into the stability problem by simply considering limiting situations. By analyzing the 
stability of these limits a qualitative picture of possible phases (or the ground state ) will emerge. In the present work 
we will investigate the effect of spin-phonon interaction in the non-adiabatic limit i.e. when phonons are certainly 
fast. In this case we will treat the phonons as fast variables and derive an effective interacting fermion model. 

We propose to study an XY spin chain whose magnetic interaction depends on the bond length. The reason for 
this study is two fold i) The undeformed spin chain Hamiltonian can be solved exactly. 

ii) The model contains the essential elements for the spin-Peierls transition, that is coupling to intermolecular 
motion. In the future it will be of interest to study the Heisenberg chain with next-nearest neighbor frustration term 
to make the model more realistic to the inorganic spin-Peierls systems. 

The paper is organized as follows. In section 2 the variational ground state energy of the XY model in presence of 
spin-phonon interaction is determined using suitable phonon states and mean field approximation taking into account 
the dimerization as well as anti-ferromagnetic ordering. The results of numerical solution and its implication are 
discussed in section 3. 

2. Formulation 

We start with the XY model in presence of spin phonon interaction on a linear chain: 
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h = j + 9(bj +k- b\ +1 - b l+1 Ms?s? +1 + srsr +1 ) +^y. b i b i (1) 

I I 
where I denotes the site index of the N site linear chain, 5 ; x and ST are components of the local X-Y spin of value 
bi(b\) is the annihilation (creation) operator for a vibration of molecule at site / and J is the magnetic exchange 
interaction between the nearest neighbor spins. Here u> accounts for the dispersionless vibrational spectra for molecular 
motion along the chain direction and g is the spin-phonon interaction [ pl| . 

We transform the spin operator to a spinless pseudo fermion representation using the Jordan Wigner transformation 

n 

l-l 

S? + iSj = S+ = exp(-i7r ]T 4^)4, (2) 

3 

l-l 

Sf~ — iSj = S[~ = exp(iir'^^ djdj)di, (3) 



Sf = l + <tfd, (4) 

to make use of the growing understanding of one dimensional Fermi system. 

After the Jordan Wigner transformation the Hamiltonian (1) can be written in terms of fermion operators d\ and 
di as 
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where 

Pi=d\d l+l +d\ +l di. (6) 

In the adiabatic approximation the spin-phonon interaction deforms the lattice to undergo the Peierls instability. 
To take into account the lattice distortion due to spin phonon coupling in our case we choose a variational wave 
function \ip} p h = U\0) with 

U = exp(±J2( b l- b M P i- p i-i)) CO 

for the phonon subsystem, where |0) is the zero phonon state and U describes a modified Lang-Firsov transformation 
p3| , ^4f . In this formalism the effective fermion Hamiltonian is written as 

H eff = (0\H T \0), (8) 

with 
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H T = U~ X HU 

= ^E + 9 -=r± ^ + b ^ P ^ 

i i 



A 9 12 4 f' 2 E nini+i + g ' 2 E^ 1 ~ 2 "0(4-i^+i + 4+i*-i) 
i i 

E^ ( 9 ) 



U) 

i 



where, 



m = d\d h (10) 

In the above equation 6/ and 6j is the creation and annihilation operator for the phonon system vibrating about the 
displaced equilibrium position ^{Pi — Pi-i) of the lattice. Clearly, A is proportional to a lattice displacement created 
by the spin-phonon interaction which has to be determined variationally. When A = g the transformation is exactly 
the Lang-Firsov p3|| transformation where the fermion-phonon term is diagonalised exactly and the fermion hopping 
term is renormalized by dressed phonons. To obtain an effective fermionic Hamiltonian we take the average over the 

A 4 

zero phonon state of the transformed phonon subsystem and neglect the terms of the order of and higher. In this 
approximation the effective Hamiltonian is 

J 3A 2 JA 2 

Heff(X) = "(I - ^) + 4+1*) + ^2 E(4*+3 + 4+3*) 

- 4g' 2 E ni + 4g' 2 E nWi+x + g' 2 E( x - 2^)(4_i*+i + 4+1 *-i) 
/ i 

+ 0(A 4 ). (12) 



Now, the above effective hamiltonian (12) is complicated enough to solve it exactly. Therefore one has to look for 
the approximate methods. We calculate the ground state energy E e ff(X) of the effective Hamiltonian (13) in the 
framework of mean field theory. We assume solutions which break the symmetry between even and odd sites with 
respect to the number of fermions on the site and to the corresponding hopping probability. Both of these kinds of 
oder open a gap at the Fermi momentum at half filling (ie. when the total spin M z is zero) and lower the ground 
state energy. We will consider four variational parameters such h e and h ( where, e implies even sites 

and o implies odd sites). All these variables are not independent because they are subject to the fermion number 
conservation constraint 

_ < 4^21 + 4 l+1 d 2l+1 > _ (n e + n ) 
n ~ 2 ~ 2 [ ' 

For half filling (M z = 0) n = i. The three remaining variational parameters are then the anti-ferromagnetic order 

parameter 



< d} 2l d 2 i - 


- dl l+1 d 2 i+i > 


(n e - n ) 
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< 4i<Wl 


- <4/+l rf 22+2 > 


(h e - h D ) 
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< 4^21+1 


+ 4/+W2/+2 > 


(h e + h a ) 
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the dimerized order parameter 

7= """"^ = v^_^ (lg) 

and average hopping probability 

Here < > implies the expectation value over the ground state. Within the limitation of Hartree-Fock approxima- 
tion our effective Hamiltonian can be written as 

H eff = |(1 - ^2) J2(d]d l+1 + d\ +1 di) + ^ ^(4^+3 + 4+3*) 



- 4. 9 ' 2 (i - m 2 - 2 7 2 )7V - 2 3 ' 2 7V + 0(4) (17) 
To diagonalise the Hamiltonian (18) we transform the operators from coordinate space to momentum space 



s f = -^E4^ as) 

c, = J=5>e-<« (19) 

Due to the reduced symmetry each fc state is coupled to the state k + n. So it is convenient to write the Hamiltonian 
in the reduced zone — ^ to ^ and label the states by band indices 1 and u. In this representation the Hamiltonian is 
a two band Hamiltonian where the band I and u are coupled to each k. 

H eff = 44 f 4 + E a i d k d i + E + E « f 4 

k k k k 
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-4 5 ' 2 (i-m 2 - 2 7 2 )7V + 0(^) (20) 
where, 

3A 2 2JA 2 

al = J(l - — T )cos(fc) + — r cos(3k) (21) 
4a; 8a; 
q\2 o /A 2 

«fc = - J (! - W fc ) - ^7j2- cos ( 3fc ) ( 22 ) 
/3 fe = -8.g' 2 (m - 2i-ysin(k)) (23) 

We diagonalise this Hamiltonian by a Bogoliubov Valatin transformation and obtain 
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k k 

- 4 5 ,2 (i - to 2 - 2 7 2 )7V (24) 



with upper and lower band 



Ef = W-f + 64--J- + WV*n(fc)» (25) 



From the equation (25) it is clear that the energy spectrum has been split into two separate bands ( for non zero m 
or 7) characterized by the Bogoliubov transformed creation (annihilation) operators a\(ak) and b\{bk)- For half-filling 
( M z = 0) the lower band is completely filled in the ground state . We will take the expectation value of the equations 
(14) (15) over the ground state to obtain a set of self-consistent equations of 

k=ir/2 



, 2 K = 7T/Z 



„ am 

fc=-7r/2 1 k 



and 



2 fc=7r/2 

g z 7 ^ sin(k) 
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k=-Ti/2 1 fel 

2 

If we retain terms up to the order of ^ the integrals involving the equations (26) and (27) can be expressed in terms 
of elliptic functions of the first kind (K(v)) and the second kind {E{v)) to give 
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with v = J g , Am . A and 

i+ .7^(1-3^ /(2^)) 



m = * v —- 2 KM (28) 

^(i-f£) 



325 '^ W) (29) 



' Trja-ii)^ 

and i/' = 1— /^(^^/^o;^)) • -^ or sman m or 7 tne equations provide asymptotic expressions for 1) the antiferromagnetic 
and 2) the dimerized phase. 

1) The antiferromagnetic phase (7 = 0; m 7^ 0) 

7 = 0; 



1.38629, 

m = pjexp{— npj{\ ) (30) 

Pj 



2) The dimerized phase (m = 0; 7 7^ 0) 



m = 0; 

pj Trpj 0.77308 , 

7= T esp(- — (1-^— ) (31) 



where pj = J ^ <- 

When m or 7 is large we solve the equation (26) or (27) numerically to get the anti-ferromagnetic or dimerized 
order parameter. 

Finally we minimize the ground state energy (Eq(X)) of the Hamiltonian -ff e // with respect to A to get A and the 
ground state energy of the system. 

3. Results and discussion 

We have computed the ground state energy, the lattice distortion and the energy gap for the XY system in pres- 
ence of spin-phonon interaction when phonon frequency is certainly fast. The results of our calculation shows some 
distinguished features in the phase diagram which is absent when phonons are treated adiabatically. In the adiabatic 
approach the spin-phonon coupling to the XY spin system always gives rise to dimerized phase where as our nona- 
diabatic approach predicts spin-liquid, dimerized and antiferromagnetic phase depending on the values of — and ^ 
ratio. Here we present the result for two values for exchange interaction — = 0.2 and 0.4. In our calculation we have 
neglected the terms of the order ^4 or higher. So we do not extend our calculation for large spin-phonon coupling 
and confine ourselves to — < 0.5 . In figure 1. we have shown the variation of lowest energy in the anti-ferromagnetic 
as well as dimerized phase to determine the ground state of the system. For — = 0.2 we find that the dimerized phase 
represents the ground state when ^ < .46 and anti-ferromagnetic phase becomes the ground state when ^ > 0.46. 
The critical value for dimerized to anti-ferromagnetic phase transition increases with the increase of ±. For J=0.4 
the above mentioned critical value of — is higher than 0.5 and it is not shown in our plot because we have confined 
our investigations within = 0.5. 

In figure 2. we show the variation of dimerized order parameter 7 with ^ for — = 0.2 and 0.4. The figure shows that 
when 2- is less than 0.1 the dimerized order parameter is almost zero. From the equation (31) we obtain 7 ~ 10 -20 
at — = 0.1. In the mean-field calculation we neglect quantum fluctuations and always find a nonzero value of the 
order parameter. However, the quantum fluctuation due to spin excitation would have drastic repercussions on the 
value of 7 and can destroy this small value of 7 to give a spin liquid phase. So one may expect a critical spin-phonon 
coupling (^j) for the on set of dimerization. Caron and Moukouri ]l4| have shown this existence of critical spin - 
phonon coupling (2s.) through the density matrix renormalization group calculation where quantum fluctuation are 
taken into account rigorously. It is also evident from our figure 2. that the dimerization order parameter decreases 
with the increase of — which means that we expect a larger ^ for — = 0.4. 

Due to Peierls instability we have two types of bond length one is greater than the unperturbed lattice constant and 
other is less than the unperturbed lattice constant. We represent this lattice deformation by < b\ + bi — bj +1 — 6; + i >= 
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(— l) l 6. In figure 3. we plot the variation of S = 2Xj with respect to £ for ^ = 0.2 and 0.4 in the dimerized phase. 
It is seen from the figure that for small £ ratio the lattice distortion is almost zero and it increases with the increase 
of - or with the decrease of — . 

We have also investigated the excitation energy gap (A) which is obtained by calculating the lowest excited state 
energy (ie. lowest energy of the upper band ) . In the figure 4. we plot the variation of A with respect to *j for 
£ = 0.2 and 0.4. Like the dimerization order parameter, A is also very small for small spin-phonon coupling. The 
result also points to the possibility of gapless spectra for £ < ^ and gapfull spectra for £ > ^. 

In figure 5. we show the variation of the critical value ^ for dimerized to antiferromagnetic phase transition with 
respect to — . We observe a decrease of — with the decrease of — and in the limit — — > the system shows an 
antiferromagnetic order for any finite spin-phonon coupling. This feature is significantly different if the spin-phonon 
coupling is treated adiabatically where the ground state of XY model with spin-lattice interaction always represents 
dimerized phase. Our results signifies that the quantum correction may play a very important role in determining the 
phases of the ground state of the spin-Peierls systems. 

In summary, we have developed a nonadiabatic approach for the interacting spin-phonon problem. We have chosen 
the Lang-Firsov variational wave function to integrate out the phonon degrees of freedom and obtained an effective 
spinless fermionic hamiltonian which is solved in the frame work of mean field approximation to calculate the ground 
state energy, the minimum excitation energy gap and the lattice distortion developed in the ground state. Our 
investigation indicates two types phase transition one is from spin liquid phase to dimerized phase and another is 
from dimerized phase to antiferromagnetic phase as we vary the spin-phonon coupling from a very low value. However, 
these phase transitions are absent if we neglect the quantum lattice fluctuations and treat the problem in the adiabatic 
approximation. So it is evident that the behavior of spin-Peierls system is significantly modified by quantum lattice 
fluctuations and an intense investigation is required to explore the effect of nonadiabatic corrections to the spin-Peierls 
systems. For a realistic calculation of inorganic spin-Peierls system like CuGeOy, one has to study the Heisenberg 
spin chain with next nearest neighbor frustration . 

Acknowledgments 

The author acknowledges the helpful discussions with E. Miiller-Hartmann, G. Bouzerar and G. Uhrig. Research was 
performed within the program of the Sonderforschungsbereich 341 supported by the Deutsche Forschungsgemeinschaft 
(DFG). 



8 



M. Hase, I. Terasaki and K. Uchinokura, Phys. Rev. Lett. 70, 3651 (1993) 

W. Palme et al. J. Appl. Phys. 79, 5384 (1996) 

J. G. Lussier et. al. J. Phys. Cond. Matter 7, L325 (1996) 

J. P. Renard et al. Europhys. Lett 30 475 (1995) . 

J. P. Boucher and L. P. Regnault J. Phys. I France 6 1, (1996). 

M. C. Cross and D. S. Fisher, Phys. Rev. B. 19 402 (1979). 

M. Fujita and K. Machida, J. Phys. Soc. J. Phys. Soc. Jpn. 53 4395 (1984). 

T. Nakano and H. Fukuyama J. Phys. Soc. Jpn. 49 1679 (1980); ibid 50 2489 (1981). 

G. S. Uhrig and H. J. Shultz, Phys. Rev. B 54 R9624 (1996); G. S. Uhrig, Phys. Rev. Lett 79 163 (1997). 
G. Castilla, S. Chakravarty, and V. J. Emery, Phys. Rev. Lett. 75 1823 (1995) 

J. Riera and A. Dorby, Phys. Rev. B. 51 16098 (1995). 

G. Bouzerar, A. P. Kamf and F. Schonfeld, preprint |cond-mat/9701 176 



G. Bouzerar, A. P. Kamf and G. I. Japaridze accepted for publication in Phys. Rev. B. ( :ond-mat/9801046) 
L. G. Caron and S. M oukouri, Phys. Rev . Lett 76 4050 (1996) 
G. S. Uhrig, preprint, |sond-mat/9801185| 



L. P. Regnault et. al., Phys. Rev. B 5579 (1996) 
M. Braden, G. Wilkendorf, J. Lorenzana, M. Ain, G. J. Mcintyre , M. Behruzi, G. Hegar, G. Dhalenne and A. Revcolvschi, 
Phys. Rev. B 54 1105 (1996). 

E. Fradkin and J. E. Hirsch Phys. Rev. B. 27 1680 (1983). 

W. P. Su, J. R. Schrieffer and A. J. Heeger Phys. Rev. B. 22 2099 (1980). 

D. K. Campbell and A. R. Bishop Nucl. Phys. B. 200 297 (1982). 

J. W. Bray et. al. Phys. Rev. Lett. 11 744 (1975). 

P. Jordan and E. Wigner, Z. Phys. 47 631 (1928). 

I. Lang and Y. A. Firsov Sov. Phys. JETP 16 1301 (1963). 

A. N. Das and S. Sil, J. Phys. Cond. Matt. 5 8265 (1993). 



9 



Figure captions 

Figure 1. Variation of the minimum energy (in the unit of the energy of the phonon) with respect to the spin- 
phonon coupling (^) in the antifcrromagnctic phase for J=0.2 (curve a) and J= 0.4 (curve c) and in the dimerized 
phase for J=0.2 (curve b) and J=0.4 (curve d). 

Figure 2. Variation of the dimerized order parameter 7 with respect to £ for J = 0.2 (curve a) and J = 0.4 (curve 
b) in the dimerized phase. 

Figure 3. Variation of the lattice distortion S with respect to £ in the dimerized phase for J = 0.2 (curve a) and 
J = 0.4 (curve b). 

Figure 4. The excitation energy gap A as a function of £ in the dimerized phase for J = 0.2 (curve a) and J = 0.4 
(curve b). 

Figure 5. Variation of the critical value ^ for dimerized to antiferromagnetic phase transition with respect to 
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